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Abstract: We present simulation results employing overlap fermions for the axial corre- 
lation functions in the e-regime of chiral perturbation theory. In this regime, finite size 
effects and topology play a dominant role. Their description by quenched chiral perturba- 
tion theory is compared to our numerical results in quenched QCD. We show that lattices 
with a linear extent L > 1.1 fm are necessary to interpret the numerical data obtained in 
distinct topological sectors in terms of the e-expansion. Such lattices are, however, still 
substantially smaller than the ones needed in standard chiral perturbation theory. How- 
ever, we also observe severe difficulties at very low values of the quark mass, in particular 
in the topologically trivial sector. 
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1. Introduction 

There are a number of questions in QCD where lattice simulations are hardly feasible in 
the physical regime, e.g., when the light quark masses assume their physical values and 
the physical volume is large enough to accommodate the light mesons. However, in some 
occasions it is possible to extract physical information from simulation results obtained 
in an unphysical setup, in particular in an unphysically small volume. The variation of 
the volume then allows for the investigation of finite size effects (FSE), which provide 
information about the same system in a large volume. Examples for this strategy are 
the non-perturbative renormalization procedure based on the Schrodinger Functional Q 
and the evaluation of a scattering length based on the analysis of the level repulsion pj. 
A further example — which will be discussed here — is the study of the dynamics of 
quasi-Goldstone bosons in the e-regime of chiral perturbation theory (xPT). 

One procedure to make use of the FSE is the transmutation of the variation of the 
volume into a variation of the physical scale. Then the effect on certain scale-dependent 
observables enables for instance the determination of a running coupling constant @. 

An alternative procedure — that we are going to deal with — is based on the fact that 
even in small volumes the FSE are characterized by infinite volume parameters. Therefore 
the evaluation of FSE provides access to physical quantities. That situation occurs in 
particular in a specific regime of xPT, where the pion Compton wave length is much larger 
than the linear extent of the volume, while the inverse linear extent is far below the QCD 
cutoff scale. The resulting FSE are described in xPT — which considers the pions as quasi- 
Goldstone bosons of the chiral symmetry breaking — by the so-called e-expansion ||]. The 
corresponding zero-modes are treated exactly by means of collective variables, while the 
non-zero modes, as well as the pion mass, are dealt with perturbatively; note that the 
fluctuations around the zero modes do fit into the box. The domain where this expansion 
applies is denoted as the e-regime. 



As a test case for QCD, the e-expansion can also be applied to the O(N) symmetric 
linear ||] and non-linear || a- model. For the case of the 4d 0(4) model the numerical tech- 
nique of extracting infinite volume parameters from FSE was applied very successfully (7j. 
In QCD such theoretical concepts are often tested in the quenched approximation where 
simulations are much easier, although the concept is designed for full QCD. It turned out 
that quenched results are quite often close to experimental quantities || , even though there 
is no reason to expect this behavior in general. At low energy, it is important that also 



quenched x?T has been worked out hence quenched QCD can be used to probe 



the applicability of this procedure. This is the motivation for our pilot study, which has 
the purpose to gain first experience on 

• the range for the physical volume where the e-expansion is applicable 

• the extent of the statistics which is required for reliable results in this framework 

• the sources of potential difficulties, and their dependence on the topological sector. 

Our statistics is rather modest and it should be enlarged for precise results. We hope, 
however, to provide a sound basis for the prospects of simulations in the e-regime, also in 
view of full QCD. 

So far the e-regime appeared as a terra incognita for numerical studies of QCD because 
until recently no formulation of lattice fermions was known which would be suitable in this 
context. The obstacles for such a formulation are linked to the problems in formulating 
chiral symmetry and topological charges v on the lattice — note that topology becomes 
extremely important for xPT in small volumes fl^| . None of these problems is solved by 
the Wilson fermion, and at small quark masses the latter also runs into technical problems 
such as the significance of exceptional configurations. The staggered fermion, on the other 
hand, seems to be "topology blind", i.e. its results do not depend on v — at least not 
for strong and moderate gauge couplings, which were tested so far — in contrast to the 
continuum |13]. 



A solution to these two problems is now available due to the recent formulation of 
chirally invariant lattice fermions. Their Dirac operator obeys the Ginsparg- Wilson rela- 



tion [15] and an operator of that kind will be used in this work. However, its application 
is very computation time demanding, hence its use in QCD is restricted to the quenched 
approximation at present and in the foreseeable future. For domain wall fermions as an- 
other realization of a Dirac operator inducing an exact lattice chiral symmetry, an attempt 
had been made to explore the e-regime [16]. However, the lattices used in that work were 



unfortunately too small for to work. 

In section |2] we describe the lattice formulation of our simulations and the formulae of 
quenched xPT which are relevant in this context. In section || our numerical results for the 
axial correlation functions in distinct topological sectors are compared to the predictions 



Using very smooth gauge fields and smeared links, topological structures can be seen by the staggered 
fermion, as a very recent study shows in the Schwinger model p4|. 



of quenched xPT. We also discuss the stability of such numerical results, which is most 
problematic in the topologically trivial sector. Section || is devoted to our conclusions 
and to an outlook on this new field of research. A synopsis of our results was given in 
refs. P'Ul- 



2. Chiral invariant lattice fermions and chiral perturbation theory 
2.1 Overlap fermions 



The lattice Dirac operator that we are going to use is the so-called overlap operator [19]. 
On a lattice of unit spacing the massless operator is given by 



D. 



(0) 



1 + A/VAtA 



D w - H , 



(2.1) 



where D\y is the Wilson-Dirac operator (with Wilson parameter 1 and k = 1/8), and fi is 
a mass parameter, which can be chosen in some interval that shrinks as the gauge coupling 
increases. Our simulations were performed at (3 = 6 and we set /i = 1.4, which is optimal 



for locality in this case 2C[. This observation referred to the Wilson gauge action that we 
were also using throughout this work. 



.Do? is a solution to the Ginsparg- Wilson relation [15] 



^75 + 75^ = jD(» l5 D(% 



(2.2) 



which excludes additive mass renormalization and exceptional configurations. Hence D, 



(0) 



has exact zero modes with a definite chirality [21]. This provides a sound definition of the 
topological charge of a lattice gauge configuration by means of the Atiyah-Singer Index 
Theorem. These amazing properties are related to the fact that Ginsparg- Wilson fermions 



have an exact lattice chiral symmetry [22], which turns into the standard chiral symmetry 
in the continuum limit. 

For practical applications we still have to insert a small bare quark mass m q . It is 
added to as follows, 

(2.3) 



The challenge to numerical studies is the inverse square root in eq. (2.1). It is basically 
equivalent — but somewhat more convenient — to evaluate the sign function of the her- 
mitean operator 75 A In our simulations we approximated this sign function by Chebyshev 
polynomials to an accuracy of 10~ 12 for the eigenvalues and 10~ 16 for the propagators. 
For the description of highly optimized algorithmic tools to deal with overlap fermions we 
refer to refs. fl2l 
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2.2 The e-regime in quenched chiral perturbation theory 

The low energy properties of QCD can be described by an effective lagrangian £ e ff, and 
^PT is a powerful tool for its construction and evaluation pjj. At leading order, C c g is 
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parametrized by the pion decay constant F n and the scalar condensate E, 



C cS [U] = ^-Tr Idf.Ud^] - l^m q Tr \ue ie/N f + U^e~ ie/N f] , (2.4) 
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where U (x) S SU(iVy) is the quasi-Goldstone boson field, Nf is the number of quark flavors 
and 9 the vacuum angle. (For simplicity we assume a single bare mass m q for all flavors.) 
Note that the values of F n and E in a finite volume and in infinite volume coincide, up to 
logarithmic corrections which are specific to the quenched approximation. 

Depending on the counting rules for the expansion, one distinguishes two regimes for 
the xPT: the "p-expansion" applies to large volumes where FSE are suppressed p6| , which 
is the standard situation. In contrast, the "e-expansion" assumes a very small quark mass 
- even below the physical masses of the u and d quarks — so that m n L <C 1, where 
m n oc y/rriq is the pion mass and L is the linear size of the the system. In this regime 
FSE are strong, but they can be described analytically for dynamical as well as quenched 
quarks J|, ||, As an important feature, the observables depend significantly on the 
topology fl~2[| , hence they should be measured by averages over gauge configurations in 
a fixed topological sector (or, more precisely, at a fixed absolute value of the topological 
charge v). In the light of the requirements of small quark masses and a well-defined 
topology, the use of Ginsparg- Wilson fermions is ideal, as we pointed out before. 

In this work we aim at a first interpretation of lattice data — obtained from simulating 
overlap fermions — in the framework of the e-expansion. In particular, our reference point 
is the formula for the axial correlation function in quenched derived in ref. [10]. 



Using this formula, the lattice data yield predictions for the leading (quenched) low energy 
constants in the effective lagrangian. In order to stay in the e-regime we have to require 
the dimensionless scaling variable 

z = m q ^V (2.5) 

to be small, i.e. clearly below 1 (where V is the volume). 2 

To render this paper self-contained we now quote the main result of quenched 
which we are going to use in our data analysis. This result was obtained |1(| consistently 
from two different methods called "supersymmetric" |27j and "replica" [28]. Ref. |l0[ 
also showed that the corresponding vector correlation function vanishes to all orders in the 
quenched e-expansion; a general argument is added in Ref. |]29| . 

We start from the definition of the bare axial current at momentum p = and Eu- 
clidean time t, 

M*) = S ^)757^(i, x) . (2.6) 

X 

The formula for the axial correlation function in a volume L 3 x T, to the first order in 



2 If we insert E = (250 MeV) 3 as a typical value from the literature, the parameters of our results in 
figures H to H correspond to z ~ 0.34. 



quenched xPT, reads 3 



(A (t) Ao(fl)) v = 2-(^ + 2m ? £ M (*)T • %(r)^ , 



hi(r) = - ( t 2 - t + - ) , t = — , 
2 V 6J T 

£ v (z) = S (z [I v (z)K v {z) + + ~) . (2.7) 

i), and K u are modified Bessel functions, z is denned in eq. (pjq), is again the topological 
charge, and /ii is a purely kinematic function. The latter shows that this correlation func- 
tion is given by a parabola in t with its minimum at T/2. This behavior is qualitatively dif- 
ferent from the infinite volume, where the correlations decay exponentially (resp. as a cosh 
function in a periodic volume). However, the small and the infinite volume have in common 
that in both cases the quenched axial correlation function to the first order only depends on 
the coupling constants F n and S. In contrast, the quenched pseudoscalar-pseudoscalar and 
the scalar-scalar correlation functions involve further parameters ||, which motivated our 



focus on the simple case of the axial correlation function (2.7). We note that the bare axial 



current in eq. (2.6) needs to be renormalized multiplicatively by applying a renormalization 



constant Za, see below. 

At this point we would like to mention that Random Matrix Theory can be applied 
to QCD also in the quenched case [|31| . Random Matrix Theory involves additional 
assumptions and it leads to explicit predictions for the distributions of the low lying eigen- 
values of the Dirac operator |3^]. If these predictions hold, they provide another way to 
extract values for £ and (somewhat more involved) for F n from lattice data. In recent 
studies of the spectrum of D$ it turned out that these Random Matrix predictions are 
essentially confirmed, if the volume is not too small, which means for the linear size: 
L > 1.2 fm fl33[ |. This agrees roughly with an earlier study using a truncated fixed point ac- 
tion on a 4 lattice [p3 . Proceeding to larger physical volumes the quality of the agreement 



with the Random Matrix prediction improves further [35]. 

In particular the probability distribution for the first non-zero Dirac eigenvalue |Ai| has 
a peak, which moves to larger values (of |Ai| VS) as \u\ increases, where v is the topological 
charge. In the neutral sector, v = 0, there is a significant density of very small eigenvalues. 



A previous study of the chiral condensate |36[] observed that this property yields severe 



problems for measurements in the sector v = 0. That observation was confirmed in a 
striking way by our investigation that we are going to present in the next section. 



3. Numerical results 



In our simulations we used the overlap operator D ov described in subsection 2.1. We 
simulated quenched QCD on periodic lattices of two sizes, 10 3 x 24 and 12 4 , both at = 6. 



3 We have an extra factor of 2 compared to ref. jlC| because our definition of the current does not involve 
the flavor space generators. 



lattice 
size 


complete 
statistics 


numbe 
v = 


r of confi 
\u\ = 1 


durations 
\u\ = 2 


10 3 x 24 


65 


20 


24 


17 


12 4 


157 


47 


78 


24 



Table 1: The statistics of our simulations on two lattices, which we simulated both at ft = 6. 




t t 

Figure 1: Data for the axial correlation function from a 10 3 x 24 lattice at j3 = 6. For comparison 
we also show the curves predicted by quenched xPT for £ = (solid line) and S = (250 MeV) 3 
(dashed line). The plot for \v\ = 1 (on the right) illustrates clearly that the lattice size chosen 
here is not compatible with ;\TT, as it could be expected also based on the comparison to Random 
Matrix Theory Q. 



This corresponds to physical volumes of (0.93 fm) 3 x 1.86 fm resp. (1.12 fm) 4 . 4 

To determine the topological charge we computed the index, i.e. the difference between 
the number of right-handed and left-handed zero modes of Dqv > by using either the Ritz 



functional method |$8| or the Arnoldi procedure |3^]. We mention in passing that the 
index was always given by the number of zero modes with only one chirality; we never 
encountered zero modes of opposite chiralities in the same configuration, as it is expected 



in general for interacting chiral fermions [24, 40 1 



Our statistics on the two lattice sizes is given in table |]. On such small volumes most 
configurations have charges \i>\ < 2, and we have the largest statistics at \u\ = 1. 

The parameter F n is determined through the minimum of the function in eq. (|2.7j ) at 
r = 1/2 and it is hence easy to fit. S, on the other hand, is related to the curvature of the 
predicted parabola and more difficult to extract, as we are going to see. 5 

First we observed also here that the physical volume must exceed the lower limit 



mentioned in subsection 2.2. As a striking example, we show that otherwise one does not 



4 Our physical scale is based on the Sommer parameter p7| . In the following we do not keep track of 
possible errors in that scale. 

5 This situation is somehow complementary to the comparison of the eigenvalue distributions to their 
Random Matrix predictions, where E is easier to determine. 
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Figure 2: Histories of the axial correlation function Ao(t)Ao(0) on a 10 3 x 24 lattice at /3 = 6 and 
uiq — 21.3 MeV. We show the histories at t = 8 and at t = 16. For some v = configurations (on 
the left) we recognize marked spikes. At \v\ = 1 (on the right) the history is significantly smoother. 
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Figure 3: The height of the two spikes of Aq{1(!))Aq{{S) in the v = history of figure g, at different 
quark masses m q = 0.2 MeV ... 21.3 MeV. 



obtain the curvature required at \v\ = 1 by any positive value of E. This is illustrated in 
figure |l| (on the right), where our data from the 10 3 x 24 lattice are shown, as well as the 
analytic curves given by eq. (|2.7l ) for S = and S = (250 MeV) 3 . The fact that those two 
curves can hardly be distinguished is typical for all sectors with v ^ 0, due to the last term 
in the expression of eq. ( |2.7| ) for 

Once we are in the e-regime, v = may look like the most convenient case from 
eq. ( |2.7D . In particular this is the only sector which does not have the problem of the tiny 
sensitivity of the predicted function on the value of S, see figure [l] on the left. Unfortunately, 
from the numerical point of view this sector is a nightmare. Its problems are related to the 
danger of very small eigenvalues, as we mentioned in subsection |2.2j . This causes the large 
error bars in figure |l| on the left. To illustrate this problem explicitly we show in Fig. || the 
histories for our measurements of ) A)(0))o and {Ao(t) Aq(0))%. In the topologically 



neutral sector there are strong spikes. We checked that these spikes appear exactly at those 
configurations which have very small eigenvalues. The height of the spikes is maximal at 
very small bare quark mass m q , as figure ^ shows. 6 In the chiral limit m q — * the spike 
height is proportional to l/|Ai|. Qualitatively this property also holds for the (smaller) 
spikes in the non-trivial sectors, where, however, the occurrence of such very small modes 
is suppressed; this suppression gets stronger for increasing values of \v\. Therefore we have 
some kind of a tuning problem for m q : on one hand it must be very small to make sure 
we are in the e-regime, but on the other hand taking it too small makes measurements 
statistically tedious, in addition to numerical problems when evaluating the propagators 
(c.f. footnote |6|). 

Returning to m q = 21.3 MeV, we would now like to get an idea of the effect of such 
spikes, i.e. we want to estimate how many configurations would be required to obtain stable 
expectation values. For this purpose we consider the contribution of the smallest (non-zero) 
eigenvalue Ai alone to the scalar condensate. This contribution reads 

^l = v[ dzlP ^ Zl) m ^ q + 2 (Zhvr ' ^ = w < m 



where we insert the probabilities P v (z\) provided by Random Matrix Theory f32fl . We 
generated fake numbers for the lowest eigenvalue with probability P u {z\) and computed 
^min f rom them. The result can only be trusted when the standard deviation becomes 
practically constant. Figure Q shows the evolution of this standard deviation as the statis- 
tics is enhanced. We infer that the sector v = requires a tremendous statistics of O(10 4 ) 
configurations, but the situation is clearly better for v ^ 0, where 0(100) configurations 
might be sufficient. 

We now present our results from the 12 4 lattice, which is close to the point where the 
applicability of should set in. The quark mass amounts again to m q = 21.3 MeV. 

Figure |H| (on the left) shows our data for the axial correlation function at \v\ = 1. 
The curve corresponding to eq. (^]7]) can be fitted well over some interval that excludes 
the points near the boundary (which are strongly affected by the contributions of excited 
states). This determines the additive constant — and therefore F n — quite accurately. 
In figure [B| on the right we plot the resulting values of F % as a function of the number 
tf of t values (around T/2) that we include in the fit. We find a decent plateau, which 
suggests a value of F n = (86.7 ± 4.0) MeV for the quenched, bare pion decay constant. 
Although the purpose of the present paper is mainly to demonstrate that measurements of 
(Ao(t)Ao(Q)) u ^lq do have the potential to determine F n , let us shortly address the question 
of renormalization. In principle, the renormalized axial correlation function 

{A*(t)A$(0)) = fi-(Ao®M°)) 

has to be considered. As an estimate of the renormalization constant Za the value de- 
termined in ref. [42] can be taken which amounts to Za — 1.55. Thus the value of the 

6 Figure ^ shows data for bare masses in the range of m, = 0.2 MeV ... 21.3 MeV. The corresponding 
overlap operators D ov could be inverted simultaneously at all these masses by means of a Multiple Mass 
Solver. However, if one involves very small masses the latter has a problem of accuracy pi]], hence we did 
not rely on it for our results at m q — 21.3 MeV and used standard inverters at this single mass only. 
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Figure 4: The standard deviation in the (fake) measurement of S^ n , the contribution of the 
lowest non-zero eigenvalue to the scalar condensate, given in eq. ( |3.l[ ). The parameters correspond 
to a 12 4 lattice at j3 = 6 and m g = 21.3 MeV, and we assumed £ = (250 MeV) 3 . 




Figure 5: The axial correlation function on a 12 4 lattice at /3 = 6 in the sector \v\ = 1 (on the left). 
The two lines — which are hardly distinguishable in this plot — are fits of the curve predicted by 
quenched %PT with the free parameter and £ = resp. (250 MeV) 3 . On the right we display 
the values of F n obtained from fits over tf points around the center at t = 6. 



pion decay constant as extracted from the bare axial correlation function receives a large 
renormalization and we find F* » 130 MeV. We stress again that the values of F n and F^ 
computed here are values of the quenched approximation. 

The difficulty to extract E from our data is illustrated in figure |5] (left). It shows fits 
to tf = 1 points with F n as a free parameter, for S fixed as resp. (250 MeV) 3 . Indeed, 
those two curves can hardly be distinguished. 

Finally, as a consistency check we also take a look at charge \v\ =2. Figure ^ shows 
the data and the predicted curve for the parameter range that we extracted before from 
\v\ = 1. We see that these curves are compatible with our data, within the errors. 
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Figure 6: Data from a 12 4 lattice at /3 = 6 in the sector |^| = 2, along with the curves from 
eq. (2.7). The dark area shows these curves for the parameters in the range that we determined 
from = 1. Wc find agreement with the \v\ = 2 data within the errors. 



4. Conclusions 

We presented a pilot study of the axial correlation function — as the simplest example 
of a meson correlation function — in the e-regime, based on quenched QCD with overlap 
fermions. We find that there are several conditions for running conclusive simulations 
with Ginsparg- Wilson fermions in the e-regime. The linear size should obey L>1.1 fm, 
and measurements have to be performed in a sector of non-trivial topology. The quark 
mass should be small for conceptual reasons, but taking it too small causes technical 
problems in the measurements. The above requirement on the linear size is consistent with 
an earlier observation based on a comparison to the Random Matrix predictions for the 
eigenvalues p3| . 

We observed that the neutral sector, v = 0, is highly unfavorable for measurements, due 
to the frequent occurrence of very small eigenvalues. The corresponding configurations show 
marked spikes, which lead to the requirement of a huge statistics of O(10 4 ) configurations. 

Therefore we worked in the sectors \u\ = 1 and 2. Although our statistics is still rather 
modest, we recognize that the data can be fitted to the functions predicted by quenched 
%PT. These fits allow for a stable determination of the quenched value of F n . With our 
present analysis we were not able to determine £ since large variations of this parameter 
hardly changes the fit to the axial correlation function. The only exception is the sector 
v = 0, see figure ||. However, exactly there the statistical results are ruined by strong 
spikes, as we pointed out before. 



As an outlook, further information could be obtained by considering in addition the 
pseudoscalar-pseudoscalar and the scalar-scalar correlation function [18]. Then the predic- 
tions involve four free parameters, hence it will be more difficult to establish a consistent 
picture. We also plan to proceed to larger physical volumes where we expect complete 
agreement with xPT. In particular, these parameters would then fix directly a one loop 
approximation to T, u in eq. (2.7), which is also a rough approximation to E. We also 
intend to verify the emerging picture with an alternative Ginsparg- Wilson fermion formu- 
lation |||. 

The results of this paper suggest that it is advantageous to work in fixed topological 
sectors with relatively large charges, say \u\ = 2... 5, in order to suppress the effects of 
low- lying (non-zero) eigenvalues. Hence algorithmic tools or modified gauge actions that 
would allow for such simulations are highly desirable and should be tested in practice to 
explore the e-regime. 
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